2000.044700/FUR 
TT3229 



Application for United States Letters Patent 

for 

RUN-TO-RUN CONTROL METHOD FOR AUTOMATED CONTROL 
OF METAL DEPOSITION PROCESSES 



m 
m 



by 

Thomas Sonderman 
Scott Bushman 
Craig William Christian 



-NUMBER:: 



m : Mkmm^ express mail receipts 



DATE OF DEPOSIT: May i 9| 2001 lllip j - 
thereby certifylthattN 

PostaliService "EXPRESS MAI l!| TO TO ADDRESS EE" servi^ nder 



37i^GJFiR. 1 :1 0 on the (dateliridicated abov^ 



addressedlfo: 




Patents, 




Ws^hingtbn^ pC 20231. 



BOX PATENT 



Signaturel 



2000.044700/FUR 
TT3229 

RUN-TO-RUN CONTROL METHOD FOR AUTOMATED CONTROL 
OF METAL DEPOSITION PROCESSES 

BACKGROUND OF THE INVENTION 

1. FIELD OF THE INVENTION 

This invention relates generally to semiconductor fabrication technology, and, more 
particularly, to a method for semiconductor fabrication supervision and optimization. 

2. DESCRIPTION OF THE RELATED ART 

There is a constant drive within the semiconductor industry to increase the quality, 
reliability and throughput of integrated circuit devices, e.g., microprocessors, memory 
devices, and the like. This drive is fueled by consumer demands for higher quality computers 
and electronic devices that operate more reliably. These demands have resulted in a continual 
improvement in the manufacture of semiconductor devices, e.g., transistors, as well as in the 
manufacture of integrated circuit devices incorporating such transistors. Additionally, 
reducing defects in the manufacture of the components of a typical transistor also lowers the 
overall cost per transistor as well as the cost of integrated circuit devices incorporating such 
transistors. 

The technologies underlying semiconductor processing tools have attracted increased 
attention over the last several years, resulting in substantial refinements. However, despite the 
advances made in this area, many of the processing tools that are currently commercially 
available suffer certain deficiencies. In particular, such tools often lack advanced process data 
monitoring capabilities, such as the ability to provide historical parametric data in a 
user-friendly format, as well as event logging, real-time graphical display of both current 
processing parameters and the processing parameters of the entire run, and remote, i.e., local 
site and worldwide, monitoring. These deficiencies can engender nonoptimal control of 
critical processing parameters, such as throughput accuracy, stability and repeatability, 
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processing temperatures, mechanical tool parameters, and the like. This variability manifests 
itself as within-run disparities, run-to-run disparities and tool-to-tool disparities that can 
propagate into deviations in product quality and performance, whereas an ideal monitoring 
and diagnostics system for such tools would provide a means of monitoring this variability, 
as well as providing means for optimizing control of critical parameters. 

Among the parameters it would be useful to monitor and control are process 
parameters related to metal deposition processing (MDP). Examples of such process 
parameters include the deposition thickness and deposition rate. The metal deposition 
processing (MDP) performance is typically affected by the consumption of the sputter target 
during consecutive processing runs. In particular, the deposition rate is affected by the sputter 
target life. Typically, the deposition rate decreases the longer the sputter target is used. 
Conventionally, this deposition rate decrease has been corrected manually by the process 
engineer, using a "seat of the pants" approach that may lead to misprocessing, e.g., the 
process may be continued for an incrementally longer duration. This may cause differences in 
wafer processing between successive runs or batches or lots of wafers, leading to decreased 
satisfactory wafer throughput, decreased reliability, decreased precision and decreased 
accuracy in the semiconductor manufacturing process. 

However, traditional statistical process control (SPC) techniques are often inadequate 
to control precisely process parameters related to metal deposition processing (MDP) in 
semiconductor and microelectronic device manufacturing so as to optimize device 
performance and yield. Typically, statistical process control (SPC) techniques set a target 
value, and a spread about the target value, for the process parameters related to metal 
deposition processing (MDP). The statistical process control (SPC) techniques then attempt 
to minimize the deviation from the target value without automatically adjusting and adapting 
the respective target values to optimize the semiconductor device performance, and/or to 
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optimize the semiconductor device yield and throughput. Furthermore, blindly minimizing 
non-adaptive processing spreads about target values may not increase processing yield and 
throughput. 

Traditional control techniques are frequently ineffective in reducing off-target 
processing and in improving sort yields. For example, the wafer electrical test (WET) 
measurements are typically not performed on processed wafers until quite a long time after 
the wafers have been processed, sometimes not until weeks later. When one or more of the 
processing steps are producing resulting wafers that wafer electrical test (WET) 
measurements indicate are unacceptable, causing the resulting wafers to be scrapped, this 
misprocessing goes undetected and uncorrected for quite a while, often for weeks, leading to 
many scrapped wafers, much wasted material and decreased overall throughput. 

The present invention is directed to overcoming, or at least reducing the effects of, 
one or more of the problems set forth above. 

SUMMARY OF THE INVENTION 

In one aspect of the present invention, a method is provided, the method comprising 
monitoring consumption of a sputter target to determine a deposition rate of a metal layer 
during metal deposition processing using the sputter target, and modeling a dependence of the 
deposition rate on at least one of deposition plasma power and deposition time. The method 
also comprises applying the deposition rate model to modify the metal deposition processing 
to form the metal layer to have a desired thickness. 

In another aspect of the present invention, a computer-readable, program storage 
device is provided, encoded with instructions that, when executed by a computer, perform a 
method, the method comprising monitoring consumption of a sputter target to determine a 
deposition rate of a metal layer during metal deposition processing using the sputter target, 
and modeling a dependence of the deposition rate on at least one of deposition plasma power 




Page 4 of 65 



2000.044700/FUR 
TT3229 

and deposition time. The method also comprises applying the deposition rate model to 
modify the metal deposition processing to form the metal layer to have a desired thickness. 

In yet another aspect of the present invention, a computer programmed to perform a 
method is provided, the method comprising monitoring consumption of a sputter target to 
determine a deposition rate of a metal layer during metal deposition processing using the 
sputter target, and modeling a dependence of the deposition rate on at least one of deposition 
plasma power and deposition time. The method also comprises applying the deposition rate 
model to modify the metal deposition processing to form the metal layer to have a desired 
thickness. 

In another aspect of the present invention, a system is provided, the system 
comprising a tool monitoring consumption of a sputter target to determine a deposition rate of 
a metal layer during metal deposition processing using the sputter target, and a computer 
modeling a dependence of the deposition rate on at least one of deposition plasma power and 
deposition time. The system also comprises a controller applying the deposition rate model to 
modify the metal deposition processing to form the metal layer to have a desired thickness. 

In yet another aspect of the present invention, a device is provided, the device 
comprising means for monitoring consumption of a sputter target to determine a deposition 
rate of a metal layer during metal deposition processing using the sputter target, and means 
for modeling a dependence of the deposition rate on at least one of deposition plasma power 
and deposition time. The device also comprises means for applying the deposition rate model 
to modify the metal deposition processing to form the metal layer to have a desired thickness. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention may be understood by reference to the following description taken in 
conjunction with the accompanying drawings, in which the leftmost significant digit(s) in the 
reference numerals denote(s) the first figure in which the respective reference numerals 
appear, and in which: 

Figures 1-19 schematically illustrate various embodiments of a method for 
manufacturing according to the present invention; and, more particularly: 

Figures 1, and 3-10 schematically illustrate a flow chart for various illustrative 
embodiments of a method according to the present invention; 

Figure 2 schematically illustrates a metal deposition rate plotted against a sputter 
target life; 

Figure 11 schematically illustrates a method for fabricating a semiconductor device 
practiced in accordance with the present invention; 

Figure 12 schematically illustrates workpieces being processed using a metal 
deposition processing (MDP) tool, using a plurality of control input signals, in accordance 
with the present invention; 

Figures 13-14 schematically illustrate one particular embodiment of the process and 
tool in Figure 12; 

Figure 15 schematically illustrates one particular embodiment of the method of 
Figure 1 1 that may be practiced with the process and tool of Figures 13-14; and 

Figures 16-19 schematically illustrate geometrically polynomial least-squares fitting, 
in accordance with the present invention. 

While the invention is susceptible to various modifications and alternative forms, 
specific embodiments thereof have been shown by way of example in the drawings and are 
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herein described in detail. It should be understood, however, that the description herein of 
specific embodiments is not intended to limit the invention to the particular forms disclosed, 
but on the contrary, the intention is to cover all modifications, equivalents, and alternatives 
falling within the spirit and scope of the invention as defined by the appended claims. 

DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS 

Illustrative embodiments of the invention are described below. In the interest of 
clarity, not all features of an actual implementation are described in this specification. It will 
of course be appreciated that in the development of any such actual embodiment, numerous 
implementation-specific decisions must be made to achieve the developers' specific goals, 
such as compliance with system-related and business-related constraints, which will vary 
from one implementation to another. Moreover, it will be appreciated that such a 
development effort might be complex and time-consuming, but would nevertheless be a 
routine undertaking for those of ordinary skill in the art having the benefit of this disclosure. 

Illustrative embodiments of a method according to the present invention are shown in 
Figures 1-19. As shown in Figure 1, a workpiece 100, such as a semiconducting substrate or 
wafer, having zero, one or more process layers and/or semiconductor devices such as a 
metal-oxide-semiconductor (MOS) transistor disposed thereon, for example, is delivered to a 
metal deposition processing (MDP) step 105. In the metal deposition processing (MDP) 
step 105, metal deposition processing may be performed on the workpiece 100. 

As shown in Figure 2, a dependence of a metal deposition rate on a sputter target life 
is illustrated. Sputter target consumption, measured by sputter target life (in arbitrary units), 
is plotted along the horizontal axis against deposition rates (in arbitrary units) plotted along 
the vertical axis. As described above, to the right of dotted line 200, generally the metal 
deposition rate 210 decreases as the sputter target life increases. The dependence of the metal 
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deposition rate on the sputter target life may be determined by modeling, as described more 
fully below. 

As shown in Figure 3, the metal deposition processing (MDP) step 105 may 
communicate with a monitoring step 110 and other processing steps 140 via bidirectional 
connections through a system communications bus 160. As shown in Figure 3, the system 
communications bus 160 also provides communications between the metal deposition 
processing (MDP) step 105, the monitoring step 110 and other processing steps 140, and an 
Advanced Process Control (APC) system 120, more fully described below. 

As shown in Figure 4, the workpiece 1 00 is sent from the metal deposition processing 
(MDP) step 105 and delivered to the monitoring step 110. In the monitoring step 1 10, one or 
more metal deposition processing (MDP) tool variables and/or one or more metal deposition 
processing (MDP) parameters during one or more metal deposition processing (MDP) runs 
may be monitored and/or measured. Examples of such tool variables and/or metal deposition 
processing (MDP) parameters may comprise the degree of sputter target consumption (as 
measured by sputter target life), deposition plasma power, deposition time, temperature, 
pressure, gas flow, other parameters that may affect (for example, increase or decrease) the 
mean-free-path of the sputtered species, and the like. As shown in Figure 4, the monitoring 
step 110 may communicate with the metal deposition processing (MDP) step 105 via the 
system communications bus 160. As shown in Figure 4, the system communications bus 160 
also provides communications between the metal deposition processing (MDP) step 105, the 
monitoring step 110, and the Advanced Process Control (APC) system 120, more fully 
described below. 

As shown in Figure 5, the workpiece 100 is sent from the monitoring step 110 and 
delivered to the other processing steps 140. In the other processing steps 140, other 
processing steps may be performed on the workpiece 100 to produce the finished 
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workpiece 100. In alternative illustrative embodiments, the workpiece 100 sent from the 
monitoring step 110 may be the finished workpiece 100, in which case, there may not be 
other processing steps 140. As shown in Figure 5, the other processing steps 140 may 
communicate with the monitoring step 110 via the system communications bus 160. As 
shown in Figure 5, the system communications bus 160 also provides communications 
between the monitoring step 110, the other processing steps 140, and the Advanced Process 
Control (APC) system 120, more fully described below. 

As shown in Figure 6, the monitored sensor data 1 1 5 is sent from the monitoring 
step 110 and delivered to the Advanced Process Control (APC) system 120. As shown in 
Figure 6, the Advanced Process Control (APC) system 120 may communicate with the 
monitoring step 110 via the system communications bus 160. Delivering the monitored 
sensor data 115 to the Advanced Process Control (APC) system 120 produces an output 
signal 125. 

As shown in Figure 7, the output signal 125 is sent from the Advanced Process 
Control (APC) system 120 and delivered to a metal deposition rate modeling with model 
inversion step 130. In the metal deposition rate modeling with model inversion step 130, the 
monitored sensor data 1 15 may be used in a metal deposition rate model, appropriate for the 
metal deposition processing (MDP) performed on the workpiece 1 00 in the metal deposition 
processing (MDP) step 105. Examples of such a metal deposition rate model, appropriate for 
metal deposition processing, may be provided by models empirically derived using MatLab®, 
Mathematica®, and the like. The use of the monitored sensor data 115 in a metal deposition 
rate model produces one or more metal deposition processing (MDP) recipe adjustments 145. 

In various illustrative embodiments, a metal deposition rate model may be built. Such 
a metal deposition rate model may also be formed by monitoring one or more metal 
deposition processing (MDP) tool variables and/or one or more metal deposition processing 
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(MDP) parameters during one or more metal deposition processing (MDP) runs. As described 
above, examples of such metal deposition processing (MDP) tool variables and/or metal 
deposition processing (MDP) parameters may comprise the degree of sputter target 
consumption (as measured by sputter target life), deposition plasma power, deposition time, 
temperature, pressure, gas flow, other parameters that may affect the mean-free-path of the 
sputtered species, and the like. In these various illustrative embodiments, building the metal 
deposition rate models may comprise fitting the collected processing data using at least one 
of polynomial curve fitting, least-squares fitting, polynomial least-squares fitting, 
non-polynomial least-squares fitting, weighted least-squares fitting, weighted polynomial 
least-squares fitting, weighted non-polynomial least-squares fitting, and the like, as described 
more fully below. 

In various illustrative embodiments, the metal deposition rate model may incorporate 
model inversion capability. For example, the dependence of the deposition rate on the 
deposition plasma power, deposition time, and/or other variables may be adequately 
represented by: deposition rate = ^{deposition plasma power, deposition time, ...). By 
inverting the model, the deposition plasma power and/or the deposition time that would be 
appropriate for a desired deposition rate may be determined. The deposition rate and the 
thickness of the metal layer may be related by: thickness = (deposition rate)(deposition time). 
Consequently, a deposited metal layer may be formed to have a desired thickness by 
inverting the deposition rate model, and choosing the deposition plasma power and/or the 
deposition time that would be appropriate for the desired deposition rate and, hence, the 
desired thickness. 

For example, the dependence of the deposition rate on the deposition plasma power 
and the deposition time may be adequately represented by: t = F(f,T) = f a T^ , where the 
deposition rate is t, the deposition plasma power is f, the deposition time is T, and a and (3 are 
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exponents, determined by fitting measured data points, as described more fully below. The 
desired thickness 0 of the metal layer may be given by: 0 -tT = F{f,T)T = f a T^ + \ for 
example. The variation in the deposition rate with variations in the deposition plasma power 
and the deposition time may be adequately represented by: St = af a ~ l T^S f + P f a T^~ l ST , 
for example. A deposited metal layer may be formed to have a desired thickness 0 by 
inverting the deposition rate model, t xlct =T^ /a f and = f^^T , and choosing the 
deposition plasma power/ (/ = 0 Ua T~ (P+1)/ot ) and/or the deposition time T 
(T = # 1/( P +I >y* -a/(p+1 >) that would be appropriate for the desired deposition rate t, and, hence, 
the desired thickness 0. 

The metal deposition rate modeling of the monitored sensor data 1 1 5 in the metal 
deposition rate modeling with model inversion step 130, may be used to alert an engineer of 
the need to adjust the processing performed in any of the processing steps, such as the metal 
deposition processing (MDP) step 105 and/or the other processing steps 140. The engineer 
may also alter and/or adjust, for example, the setpoints for the metal deposition processing 
(MDP) performed in the metal deposition processing (MDP) step 105, and/or the metal 
deposition processing (MDP) tool variable(s) and/or metal deposition processing (MDP) 
parameter(s) monitored and/or measured in the monitoring step 110. 

As shown in Figure 8, a feedback control signal 135 may be sent from the metal 
deposition rate modeling with model inversion step 130 to the metal deposition processing 
(MDP) step 105 to adjust the metal deposition processing (MDP) performed in the metal 
deposition processing (MDP) step 105. For example, a deposited metal layer may be formed 
to have a desired thickness by inverting the deposition rate model, and choosing the 
deposition plasma power and/or the deposition time that would be appropriate for the desired 
deposition rate and, hence, the desired thickness. In various alternative illustrative 
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embodiments, the feedback control signal 135 may be sent from the metal deposition rate 
modeling with model inversion step 130 to any of the other processing steps 140 to adjust the 
processing performed in any of the other processing steps 140, for example, via the system 
communications bus 160 that provides communications between the metal deposition 
processing (MDP) step 105, the monitoring step 1 10, the other processing steps 140, and the 
Advanced Process Control (APC) system 120, more fully described below. 

As shown in Figure 9, in addition to, and/or instead of, the feedback control 
signal 135, the one or more metal deposition processing (MDP) recipe adjustments 145, 
and/or an entire appropriate recipe based upon this analysis, may be sent from the metal 
deposition rate modeling with model inversion step 130 to a metal deposition processing 
(MDP) process change and control step 150. In the metal deposition processing (MDP) 
process change and control step 150, the one or more metal deposition processing (MDP) 
recipe adjustments 145 may be used in a high-level supervisory control loop. Thereafter, as 
shown in Figure 10, a feedback control signal 155 may be sent from the metal deposition 
processing (MDP) process change and control step 150 to the metal deposition processing 
(MDP) step 105 to adjust the metal deposition processing (MDP) performed in the metal 
deposition processing (MDP) step 105. For example, a deposited metal layer may be formed 
to have a desired thickness by inverting the deposition rate model, and choosing the 
deposition plasma power and/or the deposition time that would be appropriate for the desired 
deposition rate and, hence, the desired thickness. In various alternative illustrative 
embodiments, the feedback control signal 155 may be sent from the metal deposition 
processing (MDP) process change and control step 150 to any of the other processing 
steps 140 to adjust the processing performed in any of the other processing steps 140, for 
example, via the system communications bus 1 60 that provides communications between the 
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metal deposition processing (MDP) step 105, the monitoring step 110, the other processing 
steps 140, and the Advanced Process Control (APC) system 120, more fully described below. 

In various illustrative embodiments, the engineer may be provided with advanced 
process data monitoring capabilities, such as the ability to provide historical parametric data 
in a user- friendly format, as well as event logging, real-time graphical display of both current 
processing parameters and the processing parameters of the entire run, and remote, i.e., local 
site and worldwide, monitoring. These capabilities may engender more optimal control of 
critical processing parameters, such as throughput accuracy, stability and repeatability, 
processing temperatures, mechanical tool parameters, and the like. This more optimal control 
of critical processing parameters reduces this variability. This reduction in variability 
manifests itself as fewer within-run disparities, fewer run-to-run disparities and fewer 
tool-to-tool disparities. This reduction in the number of these disparities that can propagate 
means fewer deviations in product quality and performance. In such an illustrative 
embodiment of a method of manufacturing according to the present invention, a monitoring 
and diagnostics system may be provided that monitors this variability and optimizes control 
of critical parameters. 

Figure 1 1 illustrates one particular embodiment of a method 1 1 00 practiced in 
accordance with the present invention. Figure 12 illustrates one particular apparatus 1200 
with which the method 1100 may be practiced. For the sake of clarity, and to further an 
understanding of the invention, the method 1100 shall be disclosed in the context of the 
apparatus 1200. However, the invention is not so limited and admits wide variation, as is 
discussed further below. 

Referring now to both Figures 1 1 and 12, a batch or lot of workpieces or wafers 1205 
is being processed through a metal deposition processing (MDP) tool 1210. The metal 
deposition processing (MDP) tool 1210 may be any metal deposition processing (MDP) tool 
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known to the art, provided it comprises the requisite control capabilities. The metal 
deposition processing (MDP) tool 1210 comprises a metal deposition processing (MDP) tool 
controller 1215 for this purpose. The nature and function of the metal deposition processing 
(MDP) tool controller 1215 will be implementation specific. 

For instance, the metal deposition processing (MDP) tool controller 1215 may control 
metal deposition processing (MDP) control input parameters such as metal deposition 
processing (MDP) recipe control input parameters and/or setpoints. Four workpieces 1205 are 
shown in Figure 12, but the lot of workpieces or wafers, i.e., the "wafer lot," may be any 
practicable number of wafers from one to any finite number. 

The method 1100 begins, as set forth in box 1 1 20, by monitoring consumption of a 
sputter target to determine a deposition rate of a metal layer during metal deposition 
processing (MDP) performed on the workpiece 1205 in the metal deposition processing 
(MDP) tool 1210. The nature, identity, and measurement of characteristic parameters, such as 
deposition plasma power and/or deposition time and/or sputter target life, will be largely 
implementation specific and even tool specific. For instance, capabilities for monitoring 
process parameters vary, to some degree, from tool to tool. Greater sensing capabilities may 
permit wider latitude in the characteristic parameters that are identified and measured and the 
manner in which this is done. Conversely, lesser sensing capabilities may restrict this latitude. 
In turn, the metal deposition processing (MDP) control input parameters such as the metal 
deposition processing (MDP) recipe control input parameters and/or the setpoints for 
deposition plasma power and/or deposition time and/or flow rates of ambients (argon, Ar, for 
example) and/or chuck temperature and/or heat exchange temperature (for example, keeping 
one or more workpieces at approximately constant temperature, if possible) may directly 
affect the thickness of metal layers deposited on the workpiece 1205 and/or sputter target life. 
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Turning to Figure 12, in this particular embodiment, the metal deposition processing 
(MDP) process characteristic parameters are measured and/or monitored by tool sensors (not 
shown). The outputs of these tool sensors are transmitted to a computer system 1230 over a 
line 1220. The computer system 1230 analyzes these sensor outputs to identify the 
characteristic parameters. 

Returning, to Figure 1 1 , once the characteristic parameters are identified and 
measured, the method 1 100 proceeds by modeling a dependence of the deposition rate on the 
deposition plasma power and/or the deposition time using a metal deposition rate model, as 
set forth in box 1 130. As described above, examples of such a metal deposition rate model, 
appropriate for metal deposition processing, may be provided by models empirically derived 
using MatLab®, Mathematica®, and the like. The computer system 1230 in Figure 12 is, in 
this particular embodiment, programmed to model the characteristic parameter(s). The 
manner in which this modeling occurs will be implementation specific. 

In the embodiment of Figure 12, a database 1235 stores a plurality of models that 
might potentially be applied, depending upon which characteristic parameter is measured. 
This particular embodiment, therefore, requires some a priori knowledge of the characteristic 
parameters that might be measured. The computer system 1230 then extracts an appropriate 
model from the database 1235 of potential models to apply to the measured characteristic 
parameters. If the database 1235 does not include an appropriate model, then the 
characteristic parameter may be ignored, or the computer system 1230 may attempt to 
develop one, if so programmed. The database 1235 may be stored on any kind of 
computer-readable, program storage medium, such as an optical disk 1240, a floppy 
disk 1245, or a hard disk drive (not shown) of the computer system 1230. The database 1235 
may also be stored on a separate computer system (not shown) that interfaces with the 
computer system 1230. 
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Modeling of the measured characteristic parameter may be implemented differently in 
alternative embodiments. For instance, the computer system 1230 may be programmed using 
some form of artificial intelligence to analyze the sensor outputs and controller inputs to 
develop a model on-the-fly in a real-time implementation. This approach might be a useful 
adjunct to the embodiment illustrated in Figure 12, and discussed above, where characteristic 
parameters are measured and identified for which the database 1235 has no appropriate 
model. 

The method 1 100 of Figure 1 1 then proceeds by applying the deposition rate model to 
modify the metal deposition processing to form the metal layer to have a desired thickness, as 
set forth in box 1 140. Depending on the implementation, applying the model may yield either 
a new value for a metal deposition processing (MDP) control input parameter or a correction 
to an existing metal deposition processing (MDP) control input parameter. In various 
illustrative embodiments, a multiplicity of control input recipes may be stored and an 
appropriate one of these may be selected based upon one or more of the determined 
parameters. The new metal deposition processing (MDP) control input is then formulated 
from the value yielded by the model and is transmitted to the metal deposition processing 
(MDP) tool controller 1215 over the line 1220. The metal deposition processing (MDP) tool 
controller 1215 then controls subsequent metal deposition processing (MDP) process 
operations in accordance with the new metal deposition processing (MDP) control inputs. 

Some alternative embodiments may employ a form of feedback to improve the 
modeling of characteristic parameters. The implementation of this feedback is dependent on 
several disparate facts, including the tool's sensing capabilities and economics. One 
technique for doing this would be to monitor at least one effect of the model's 
implementation and update the model based on the effect(s) monitored. The update may also 
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depend on the model. For instance, a linear model may require a different update than would 
a non-linear model, all other factors being the same. 

As is evident from the discussion above, some features of the present invention may 
be implemented in software. For instance, the acts set forth in the boxes 1120-1140 in 
Figure 1 1 are, in the illustrated embodiment, software-implemented, in whole or in part. 
Thus, some features of the present invention are implemented as instructions encoded on a 
computer-readable, program storage medium. The program storage medium may be of any 
type suitable to the particular implementation. However, the program storage medium will 
typically be magnetic, such as the floppy disk 1245 or the computer 1230 hard disk drive (not 
shown), or optical, such as the optical disk 1240. When these instructions are executed by a 
computer, they perform the disclosed functions. The computer may be a desktop computer, 
such as the computer 1230. However, the computer might alternatively be a processor 
embedded in the metal deposition processing (MDP) tool 1210. The computer might also be a 
laptop, a workstation, or a mainframe in various other embodiments. The scope of the 
invention is not limited by the type or nature of the program storage medium or computer 
with which embodiments of the invention might be implemented. 

Thus, some portions of the detailed descriptions herein are, or may be, presented in 
terms of algorithms, functions, techniques, and/or processes. These terms enable those skilled 
in the art most effectively to convey the substance of their work to others skilled in the art. 
These terms are here, and are generally, conceived to be a self-consistent sequence of steps 
leading to a desired result. The steps are those requiring physical manipulations of physical 
quantities. Usually, though not necessarily, these quantities take the form of electromagnetic 
signals capable of being stored, transferred, combined, compared, and otherwise manipulated. 

It has proven convenient at times, principally for reasons of common usage, to refer to 
these signals as bits, values, elements, symbols, characters, terms, numbers, and the like. All 
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of these and similar terms are to be associated with the appropriate physical quantities and are 
merely convenient labels applied to these quantities and actions. Unless specifically stated 
otherwise, or as may be apparent from the discussion, terms such as "processing," 
"computing," "calculating," "determining," "displaying," and the like, used herein refer to the 
action(s) and processes of a computer system, or similar electronic and/or mechanical 
computing device, that manipulates and transforms data, represented as physical 
(electromagnetic) quantities within the computer system's registers and/or memories, into 
other data similarly represented as physical quantities within the computer system's 
memories and/or registers and/or other such information storage, transmission and/or display 
devices. 

Construction of an Illustrative Apparatus. An exemplary embodiment 1300 of the 
apparatus 1200 in Figure 12 is illustrated in Figures 13-14, in which the apparatus 1300 
comprises a portion of an Advanced Process Control ("APC") system. Figures 13-14 are 
conceptualized, structural and functional block diagrams, respectively, of the apparatus 1300. 
A set of processing steps is performed on a lot of workpieces 1305 on a metal deposition 
processing (MDP) tool 1310. Because the apparatus 1300 is part of an Advanced Process 
Control (APC) system, the workpieces 1305 are processed on a run-to-run basis. Thus, 
process adjustments are made and held constant for the duration of a run, based on run-level 
measurements or averages. A "run" may be a lot, a batch of lots, or even an individual wafer. 

In this particular embodiment, the workpieces 1305 are processed by the metal 

deposition processing (MDP) tool 1310 and various operations in the process are controlled 

by a plurality of metal deposition processing (MDP) control input signals on a line 1320 

between the metal deposition processing (MDP) tool 1310 and a workstation 1330. 

a 

Exemplary metal deposition processing (MDP) control inputs for this embodiment might 
include those for the setpoints for deposition plasma power, deposition time, temperature, 
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pressure, gas flow, other parameters that may affect (for example, increase or decrease) the 
mean-free-path of the sputtered species, and the like. 

When a process step in the metal deposition processing (MDP) tool 1310 is 
concluded, the semiconductor workpieces 1305 being processed in the metal deposition 
processing (MDP) tool 1310 are examined at a review station 1317. The metal deposition 
processing (MDP) control inputs generally affect the characteristic parameters of the 
semiconductor workpieces 1305 measured at the review station 1317, and, hence, the 
variability and properties of the acts performed by the metal deposition processing (MDP) 
tool 1310 on the workpieces 1305. Once errors are determined from the examination after the 
run of a lot of workpieces 1305, the metal deposition processing (MDP) control inputs on the 
line 1320 are modified for a subsequent run of a lot of workpieces 1305. Modifying the 
control signals on the line 1320 is designed to improve the next processing performed by the 
metal deposition processing (MDP) tool 1310. The modification is performed in accordance 
with one particular embodiment of the method 1 1 00 set forth in Figure 1 1 , as described more 
fully below. Once the relevant metal deposition processing (MDP) control input signals for 
the metal deposition processing (MDP) tool 1310 are updated, the metal deposition 
processing (MDP) control input signals with new settings are used for a. subsequent run of 
semiconductor devices. 

Referring now to both Figures 1 3 and 1 4, the metal deposition processing (MDP) 
tool 1310 communicates with a manufacturing framework comprising a network of 
processing modules. One such module is an Advanced Process Control (APC) system 
manager 1440 resident on the computer 1340. This network of processing modules 
constitutes the Advanced Process Control (APC) system. The metal deposition processing 
(MDP) tool 1310 generally comprises an equipment interface 1410 and a sensor 
interface 1415. A machine interface 1430 resides on the workstation 1330. The machine 
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interface 1430 bridges the gap between the Advanced Process Control (APC) framework, 
e.g., the Advanced Process Control (APC) system manager 1440, and the equipment 
interface 1410. Thus, the machine interface 1430 interfaces the metal deposition processing 
(MDP) tool 1310 with the Advanced Process Control (APC) framework and supports 
machine setup, activation, monitoring, and data collection. The sensor interface 1415 
provides the appropriate interface environment to communicate with external sensors such as 
LabView® or other sensor bus-based data acquisition software. Both the machine 
interface 1430 and the sensor interface 1415 use a set of functionalities (such as a 
communication standard) to collect data to be used. The equipment interface 1410 and the 
sensor interface 1415 communicate over the line 1320 with the machine interface 1430 
resident on the workstation 1330. 

More particularly, the machine interface 1430 receives commands, status events, and 
collected data from the equipment interface 1410 and forwards these as needed to other 
Advanced Process Control (APC) components and event channels. In turn, responses from 
Advanced Process Control (APC) components are received by the machine interface 1430 
and rerouted to the equipment interface 1410. The machine interface 1430 also reformats and 
restructures messages and data as necessary. The machine interface 1430 supports the 
startup/shutdown procedures within the Advanced Process Control (APC) System 
Manager 1440. It also serves as an Advanced Process Control (APC) data collector, buffering 
data collected by the equipment interface 1410, and emitting appropriate data collection 
signals. 

In the particular embodiment illustrated, the Advanced Process Control (APC) system 
is a factory-wide software system, but this is not necessary to the practice of the invention. 
The control strategies taught by the present invention can be applied to virtually any 
semiconductor metal deposition processing (MDP) tool on a factory floor. Indeed, the present 
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invention may be simultaneously employed on multiple metal deposition processing (MDP) 
tools in the same factory or in the same fabrication process. The Advanced Process Control 
(APC) framework permits remote access and monitoring of the process performance. 
Furthermore, by utilizing the Advanced Process Control (APC) framework, data storage can 
be more convenient, more flexible, and less expensive than data storage on local drives. 
However, the present invention may be employed, in some alternative embodiments, on local 
drives. 

The illustrated embodiment deploys the present invention onto the Advanced Process 
Control (APC) framework utilizing a number of software components. In addition to 
components within the Advanced Process Control (APC) framework, a computer script is 
written for each of the semiconductor metal deposition processing (MDP) tools involved in 
the control system. When a semiconductor metal deposition processing (MDP) tool in the 
control system is started in the semiconductor manufacturing fab, the semiconductor metal 
deposition processing (MDP) tool generally calls upon a script to initiate the action that is 
required by the metal deposition processing (MDP) tool controller. The control methods are 
generally defined and performed using these scripts. The development of these scripts can 
comprise a significant portion of the development of a control system. 

In this particular embodiment, there are several separate software scripts that perform 
the tasks involved in controlling the metal deposition processing (MDP) operation. There is 
one script for the metal deposition processing (MDP) tool 1310, including the review 
station 1317 and the metal deposition processing (MDP) tool controller 1315. There is also a 
script to handle the actual data capture from the review station 1317 and another script that 
contains common procedures that can be referenced by any of the other scripts. There is also 
a script for the Advanced Process Control (APC) system manager 1440. The precise number 
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of scripts, however, is implementation specific and alternative embodiments may use other 
numbers of scripts. 

Operation of an Illustrative Apparatus. Figure 15 illustrates one particular 
embodiment 1500 of the method 1 100 in Figure 11. The method 1500 may be practiced with 
the apparatus 1300 illustrated in Figures 13-14, but the invention is not so limited. The 
method 1500 may be practiced with any apparatus that may perform the functions set forth in 
Figure 15. Furthermore, the method 1100 in Figure 11 may be practiced in embodiments 
alternative to the method 1500 in Figure 15. 

Referring now to all of Figures 13-15, the method 1500 begins with processing a lot 
of workpieces 1305 through a metal deposition processing (MDP) tool, such as the metal 
deposition processing (MDP) tool 1310, as set forth in box 1510. In this particular 
embodiment, the metal deposition processing (MDP) tool 1310 has been initialized for 
processing by the Advanced Process Control (APC) system manager 1440 through the 
machine interface 1430 and the equipment interface 1410. In this particular embodiment, 
before the metal deposition processing (MDP) tool 1310 is run, the Advanced Process 
Control (APC) system manager script is called to initialize the metal deposition processing 
(MDP) tool 1310. At this step, the script records the identification number of the metal 
deposition processing (MDP) tool 1310 and the lot number of the workpieces 1305. The 
identification number is then stored against the lot number in a data store 1360. The rest of 
the script, such as the APCData call and the Setup and StartMachine calls, are formulated 
with blank or dummy data in order to force the machine to use default settings. 

As part of this initialization, the initial setpoints for metal deposition processing 
(MDP) control are provided to the metal deposition processing (MDP) tool controller 1315 
over the line 1320. These initial setpoints may be determined and implemented in any 
suitable manner known to the art. In the particular embodiment illustrated, metal deposition 
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processing (MDP) controls are implemented by control threads. Each control thread acts like 
a separate controller and is differentiated by various process conditions. For metal deposition 
processing (MDP) control, the control threads are separated by a combination of different 
conditions. These conditions may include, for example, the semiconductor metal deposition 
processing (MDP) tool 1310 currently processing the wafer lot, the semiconductor product, 
the semiconductor manufacturing operation, and one or more of the semiconductor 
processing tools (not shown) that previously processed the semiconductor wafer lot. 

Control threads are separated because different process conditions affect the metal 
deposition processing (MDP) error differently. By isolating each of the process conditions 
into its own corresponding control thread, the metal deposition processing (MDP) error can 
become a more accurate portrayal of the conditions in which a subsequent semiconductor 
wafer lot in the control thread will be processed. Since the error measurement is more 
relevant, changes to the metal deposition processing (MDP) control input signals based upon 
the error will be more appropriate. 

The control thread for the metal deposition processing (MDP) control scheme 
depends upon the current metal deposition processing (MDP) tool, current operation, the 
product code for the current lot, and the identification number at a previous processing step. 
The first three parameters are generally found in the context information that is passed to the 
script from the metal deposition processing (MDP) tool 1310. The fourth parameter is 
generally stored when the lot is previously processed. Once all four parameters are defined, 
they are combined to form the control thread name; MDPP02_OPER01_PROD01_MDPP01 
is an example of a control thread name. The control thread name is also stored in 
correspondence to the wafer lot number in the data store 1360. 

Once the lot is associated with a control thread name, the initial settings for that 
control thread are generally retrieved from the data store 1360. There are at least two 
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possibilities when the call is made for the information. One possibility is that there are no 
settings stored under the current control thread name. This can happen when the control 
thread is new, or if the information was lost or deleted. In these cases, the script initializes the 
control thread assuming that there is no error associated with it and uses the target values of 
the metal deposition processing (MDP) as the metal deposition processing (MDP) control 
input settings. It is preferred that the controllers use the default machine settings as the initial 
settings. By assuming some settings, metal deposition processing (MDP) errors can be related 
back to the control settings to facilitate feedback control. 

Another possibility is that the initial settings are stored under the control thread name. 
In this case, one or more wafer lots have been processed under the same control thread name 
as the current wafer lot, and have also been measured for metal deposition processing (MDP) 
error using the review station 1317. When this information exists, the metal deposition 
processing (MDP) control input signal settings are retrieved from the data store 1360. These 
settings are then downloaded to the metal deposition processing (MDP) tool 1310. 

The workpieces 1305 are processed through the metal deposition processing (MDP) 
tool 1310. This comprises, in the embodiment illustrated, subjecting the workpieces 1305 to a 
metal deposition process. The workpieces 1305 are measured on the review station 1317 after 
their metal deposition processing (MDP) on the metal deposition processing (MDP) 
tool 1310. The review station 1317 examines the workpieces 1305 after they are processed 
for a number of errors. The data generated by the instruments of the review station 1317 is 
passed to the machine interface 1430 via sensor interface 1415 and the line 1320. The review 
station script begins with a number of Advanced Process Control (APC) commands for the 
collection of data. The review station script then locks itself in place and activates a data 
available script. This script facilitates the actual transfer of the data from the review 
station 1317 to the Advanced Process Control (APC) framework. Once the transfer is 
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completed, the script exits and unlocks the review station script. The interaction with the 
review station 1317 is then generally complete. 

As will be appreciated by those skilled in the art having the benefit of this disclosure, 
the data generated by the review station 1317 should be preprocessed for use. Review 
stations, such as KLA review stations, provide the control algorithms for measuring the 
control error. Each of the error measurements, in this particular embodiment, corresponds to 
one of the metal deposition processing (MDP) control input signals on the line 1320 in a 
direct manner. Before the error can be utilized to correct the metal deposition processing 
(MDP) control input signal, a certain amount of preprocessing is generally completed. 

For example, preprocessing may include outlier rejection. Outlier rejection is a gross 
error check ensuring that the received data is reasonable in light of the historical performance 
of the process. This procedure involves comparing each of the metal deposition processing 
(MDP) errors to its corresponding predetermined boundary parameter. In one embodiment, 
even if one of the predetermined boundaries is exceeded, the error data from the entire 
semiconductor wafer lot is generally rejected. 

To determine the limits of the outlier rejection, thousands of actual semiconductor 
manufacturing fabrication ("fab") data points are collected. The standard deviation for each 
error parameter in this collection of data is then calculated. In one embodiment, for outlier 
rejection, nine times the standard deviation (both positive and negative) is generally chosen 
as the predetermined boundary. This was done primarily to ensure that only the points that 
are significantly outside the normal operating conditions of the process are rejected. 

Preprocessing may also smooth the data, which is also known as filtering. Filtering is 
important because the error measurements are subject to a certain amount of randomness, 
such that the error significantly deviates in value. Filtering the review station data results in a 
more accurate assessment of the error in the metal deposition processing (MDP) control input 
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signal settings. In one embodiment, the metal deposition processing (MDP) control scheme 
utilizes a filtering procedure known as an Exponentially-Weighted Moving Average 
("EWMA") filter, although other filtering procedures can be utilized in this context. 
One embodiment for the EWMA filter is represented by Equation (1): 
AVG N = W * M c + (1-W) * AVGp (1) 
where 

AVGn = the new EWMA average; 
W = a weight for the new average (AVG N ); 
Mc = the current measurement; and 
AVGp = the previous EWMA average. 

The weight is an adjustable parameter that can be used to control the amount of 
filtering and is generally between zero and one. The weight represents the confidence in the 
accuracy of the current data point. If the measurement is considered accurate, the weight 
should be close to one. If there were a significant amount of fluctuations in the process, then 
a number closer to zero would be appropriate. 

In one embodiment, there are at least two techniques for utilizing the EWMA filtering 
process. The first technique uses the previous average, the weight, and the current 
measurement as described above. Among the advantages of utilizing the first implementation 
are ease of use and minimal data storage. One of the disadvantages of utilizing the first 
implementation is that this method generally does not retain much process information. 
Furthermore, the previous average calculated in this manner would be made up of every data 
point that preceded it, which may be undesirable. The second technique retains only some of 
the data and calculates the average from the raw data each time. 

The manufacturing environment in the semiconductor manufacturing fab presents 
some unique challenges. The order that the semiconductor wafer lots are processed through 
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an metal deposition processing (MDP) tool may not correspond to the order in which they are 
read on the review station. This could lead to the data points being added to the EWMA 
average out of sequence. Semiconductor wafer lots may be analyzed more than once to verify 
the error measurements. With no data retention, both readings would contribute to the 
EWMA average, which may be an undesirable characteristic. Furthermore, some of the 
control threads may have low volume, which may cause the previous average to be outdated 
such that it may not be able to accurately represent the error in the metal deposition 
processing (MDP) control input signal settings. 

The metal deposition processing (MDP) tool controller 1315, in this particular 
embodiment, uses limited storage of data to calculate the EWMA filtered error, i.e., the first 
technique. Wafer lot data, including the lot number, the time the lot was processed, and the 
multiple error estimates, are stored in the data store 1360 under the control thread name. 
When a new set of data is collected, the stack of data is retrieved from data store 1360 and 
analyzed. The lot number of the current lot being processed is compared to those in the stack. 
If the lot number matches any of the data present there, the error measurements are replaced. 
Otherwise, the data point is added to the current stack in chronological order, according to the 
time periods when the lots were processed. In one embodiment, any data point within the 
stack that is over 128 hours old is removed. Once the aforementioned steps are complete, the 
new filter average is calculated and stored to data store 1360. 

Thus, the data is collected and preprocessed, and then processed to generate an 
estimate of the current errors in the metal deposition processing (MDP) control input signal 
settings. First, the data is passed to a compiled MatLab® plug-in that performs the outlier 
rejection criteria described above. The inputs to a plug-in interface are the multiple error 
measurements and an array containing boundary values. The return from the plug-in interface 
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is a single toggle variable. A nonzero return denotes that it has failed the rejection criteria, 
otherwise the variable returns the default value of zero and the script continues to process. 

After the outlier rejection is completed, the data is passed to the EWMA filtering 
procedure. The controller data for the control thread name associated with the lot is retrieved, 
and all of the relevant operation upon the stack of lot data is carried out. This comprises 
replacing redundant data or removing older data. Once the data stack is adequately prepared, 
it is parsed into ascending time-ordered arrays that correspond to the error values. These 
arrays are fed into the EWMA plug-in along with an array of the parameter required for its 
execution. In one embodiment, the return from the plug-in is comprised of the six filtered 
error values. 

Returning to Figure 15, data preprocessing comprises monitoring the target life of a 
sputter target to determine a deposition rate of a metal layer during metal deposition 
processing of the workpiece 1305 in the metal deposition processing (MDP) tool 1310 
variables, as set forth in box 1520. Known, potential characteristic parameters, such as the 
target life of a sputter target, may be identified by characteristic data patterns or may be 
identified as known consequences of modifications to metal deposition processing (MDP) 
control. In turn, the metal deposition processing (MDP) control input parameters such as the 
metal deposition processing (MDP) recipe control input parameters and/or the setpoints for 
deposition plasma power and/or deposition time and/or flow rates of ambients and/or chuck 
temperature and/or heat exchange temperature may directly affect the sputter target life 
and/or the thickness of metal layers deposited on the workpiece 1205. 

The next step in the control process is to calculate the new settings for the metal 
deposition processing (MDP) tool controller 1315 of the metal deposition processing (MDP) 
tool 1310. The previous settings for the control thread corresponding to the current wafer lot 
are retrieved from the data store 1360. This data is paired along with the current set of metal 
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deposition processing (MDP) errors. The new settings are calculated by calling a compiled 
MatLab® plug-in. This application incorporates a number of inputs, performs calculations in 
a separate execution component, and returns a number of outputs to the main script. 
Generally, the inputs of the MatLab® plug-in are the metal deposition processing (MDP) 
control input signal settings, the review station 1317 errors, an array of parameters that are 
necessary for the control algorithm, and a currently unused flag error. The outputs of the 
MatLab® plug-in are the new controller settings, calculated in the plug-in according to the 
controller algorithm described above. 

A metal deposition processing (MDP) process engineer or a control engineer, who 
generally determines the actual form and extent of the control action, can set the parameters. 
They include the threshold values, maximum step sizes, controller weights, and target values. 
Once the new parameter settings are calculated, the script stores the setting in the data 
store 1360 such that the metal deposition processing (MDP) tool 1310 can retrieve them for 
the next wafer lot to be processed. The principles taught by the present invention can be 
implemented into other types of manufacturing frameworks. 

Returning again to Figure 15, the calculation of new settings comprises, as set forth in 
box 1530, modeling the dependence of the deposition rate on the deposition plasma power 
and/or the deposition time characteristic parameter(s) using a metal deposition rate model. 
This modeling may be performed by the MatLab® plug-in. In this particular embodiment, 
only known, potential characteristic parameters are modeled and the models are stored in a 
database 1335 accessed by a machine interface 1430. The database 1335 may reside on the 
workstation 1330, as shown, or some other part of the Advanced Process Control (APC) 
framework. For instance, the models might be stored in the data store 1360 managed by the 
Advanced Process Control (APC) system manager 1440 in alternative embodiments. The 
model will generally be a mathematical model, i.e., an equation describing how the change(s) 
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in metal deposition processing (MDP) recipe control(s) affects the metal deposition 
processing (MDP) performance, and the like. The models described in various illustrative 
embodiments given above, and described more fully below, are examples of such models. As 
described above, examples of such metal deposition rate models, appropriate for metal 
deposition processing, may be provided by models empirically derived using MatLab®, 
Mathematical and the like. 

The particular model used will be implementation specific, depending upon the 
particular metal deposition processing (MDP) tool 1310 and the particular characteristic 
parameter(s) being modeled. Whether the relationship in the model is linear or non-linear will 
be dependent on the particular parameter(s) involved. 

The new settings are then transmitted to and applied by the metal deposition 
processing (MDP) tool controller 1315. Thus, returning now to Figure 15, once the 
characteristic parameter(s) are modeled, the deposition rate model is applied to modify the 
metal deposition processing recipe parameter(s) to form the metal layer to have a desired 
thickness, described more fully above, as set forth in box 1540. In this particular 
embodiment, the machine interface 1430 retrieves the model from the database 1335, plugs in 
the respective value(s), and determines the necessary change(s) in the metal deposition 
processing (MDP) recipe control input parameter(s). The change is then communicated by 
the machine interface 1430 to the equipment interface 1410 over the line 1320. The 
equipment interface 1410 then implements the change. 

The present embodiment furthermore provides that the models be updated. This 
comprises, as set forth in boxes 1550-1560 of Figure 15, monitoring the effect of modifying 
the metal deposition processing (MDP) recipe control input parameter(s) (box 1550) and 
updating the applied model(s) (box 1560) based on the effect(s) monitored. For instance, 
various aspects of the operation of the metal deposition processing (MDP) tool 1310 will 
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change as the metal deposition processing (MDP) tool 1310 ages. By monitoring the effect of 
the metal deposition processing (MDP) recipe change(s) implemented as a result of the 
characteristic parameter measurement, the value(s) could be updated to yield improved 
performance. 

As noted above, this particular embodiment implements an Advanced Process Control 
(APC) system. Thus, changes are implemented "between" lots. The actions set forth in the 
boxes 1520-1560 are implemented after the current lot is processed and before the second lot 
is processed, as set forth in box 1570 of Figure 15. However, the invention is not so limited. 
Furthermore, as noted above, a lot may constitute any practicable number of wafers from one 
to several thousand (or practically any finite number). What constitutes a "lot" is 
implementation specific, and so the point of the fabrication process in which the updates 
occur will vary from implementation to implementation. 

In various illustrative embodiments, the dependence of the deposition rate on the 
target life of the sputter target, the deposition plasma power, the deposition time, and the like, 
may be determined by modeling and/or fitting previously obtained metal deposition 
processing data. In these various illustrative embodiments, building the models may comprise 
fitting the collected processing data using at least one of polynomial curve fitting, 
least-squares fitting, polynomial least-squares fitting, non-polynomial least-squares fitting, 
weighted least-squares fitting, weighted polynomial least-squares fitting, and weighted 
non-polynomial least-squares fitting, and the like. 

In various illustrative embodiments, samples may be collected for 7V+1 data points 
(*/0>/)> where / = 1, 2, N 9 JV+1, and a polynomial of degree N 9 

N 

P N (x) = a 0 + a x x + a 2 x 2 +-- + a k x k H ha N x N =^a k x k , may be fit to the 7V+1 data points 

(Xj,yi). For example, 100 time data points (N= 99) may be taken relating the degree of sputter 
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target consumption, as measured by sputter target life p, the deposition plasma power f, 
and/or the deposition time 7, during one or more metal deposition processing steps, to the 
deposition rate t of metal layers being formed during the respective one or more metal 
deposition processing steps, resulting in respective sets of N+l data points (piJi), (/J,f,-), and/or 
(TiJi). The values may be the actually measured values of the metal deposition processing 
tool variables and/or metal deposition processing parameters, or ratios of actually measured 
values (normalized to respective reference setpoints), or logarithms of such ratios, for 
example. Polynomial interpolation is described, for example, in Numerical Methods for 
Scientists and Engineers, by R.W. Hamming, Dover Publications, New York, 1986, at pages 
230-235. The requirement that the polynomial Pn (x) pass through the N+l data points (x h yi) 

N 

is y i = P N (Xf) = ^a k X; , for z = 1, 2, N 9 N+l, a set of N+l conditions. These N+l 

*=o 

conditions then completely determine the N+l coefficients for k=0, 1, . . ., N. 

The determinant of the coefficients of the unknown coefficients au is the 

|x*|, where /= 1, 2, N+l, 

and k=0, 1, . N. The Vandermonde determinant V^+\, considered as a function of the 
variables x h Vn+i = Vn+i(x\jc2, - • • r^NjCN+i), is clearly a polynomial in the variables x„ as may 
be seen by expanding out the determinant, and a count of the exponents shows that the degree 

* N(N + l) 

of the polynomial is 0 + 1 + 2 + 3 + *•* + k H — + N = 2^k = — (for example, the 

diagonal term of the Vandermonde determinant Vn+i is 1 ■ x 2 -x\ x^" 1 • x^ +1 ). 



Vandermonde determinant: V N+l = 



X 
X 



N 
N 
N 

N+l 
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Now, if jc/v+i =xj 9 for 7 = 1, 2, N 9 then the Vandermonde determinant V N +\ = 0, 
since any determinant with two identical rows vanishes, so the Vandermonde determinant 
V N +\ must have the factors (x N +\-xj), for 7=1, 2, N 9 corresponding to the yV factors 

n(*jv+i ~~ x j)' Similarly, if x N = x j9 for j = 1, 2, N-l 9 then the Vandermonde determinant 
Vn+\ = 0, so the Vandermonde determinant Vn+\ must also have the factors (xyv-x 7 ), for 7 = 1,2, 

JV-1 

. . . , yV-1, corresponding to the N-l factors Y\ \ x n ~ x j)- Generally, if x m = x j9 for 7 < m 9 where 

/w = 2 ? /V, vV+1, then the Vandermonde determinant V N +\ = 0, so the Vandermonde 
determinant F/v+i must have all the factors (x m -xj) 9 for J<m, where m = 2 9 N 9 N+\ 9 

corresponding to the factors Y\ \ x m ~ x j)- Altogether, this represents a polynomial of degree 

m> 7=1 

A N(N + Y) 

/V + (vV-l) + »- +fc + -- - + 2 + l = ^A: = — , since, when m = /V+l, for example,/ may 

*=i 2 

take on any of /V values, 7 = 1 , 2, . . N 9 and when m = N 9 j may take on any of yV-1 values, 
7 = 1, 2, . . N-l 9 and so forth (for example, when m = 3 9 j may take only two values, 7 = 1,2, 
and when = 2,7 may take only one value,/ = 1), which means that all the factors have been 
accounted for and all that remains is to find any multiplicative constant by which these two 
representations for the Vandermonde determinant V N +\ might differ. As noted above, the 

diagonal term of the Vandermonde determinant V N +\ is 1 • x 2 • x] x* -1 • x^ +1 , and this may 

be compared to the term from the left-hand sides of the product of factors 

Ilk. -*;)=n( X " + l - X j)Tl( X N ~ X j)'"Tl( X l - X j)F[{ X 2 ~ X j)> X N+X ' X N~ l A '*2> 
m>j=\ 7=1 7=1 7=1 7=1 

which is identical, so the multiplicative constant is unity and the Vandermonde determinant 



Vn+i is V N+l (x l9 x 29 '- 9 x N+l )=\xf\ = Y\{x m -Xj) 



m>j=\ 
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This factorization of the Vandermonde determinant V N+ \ shows that if x, * x j9 for / ^ j, 
then the Vandermonde determinant V N +\ cannot be zero, which means that it is always 
possible to solve for the unknown coefficients since the Vandermonde determinant V^+x is 
the determinant of the coefficients of the unknown coefficients Solving for the unknown 
coefficients using determinants, for example, substituting the results into the polynomial 

N 

of degree N 9 P N (x) = ^a k x k , and rearranging suitably gives the determinant equation 

k=0 



y i 
y> i 
y 2 i 


x l 
x 2 


x 2 ■ 
x 2 ■ 
x 2 ■ 


■■ x N 
■ X? 

" 


- 0 , which is the solution to the polynomial fit. This may be 


y N 1 
y N+ i 1 


x N 

X N+l 


x 2 ■ 

N 

x 2 


■• x N 





seen directly as follows. Expanding this determinant by the elements of the top row, this is 
clearly a polynomial of degree N. The coefficient of the element y in the first row in the 
expansion of this determinant by the elements of the top row is none other than the 
Vandermonde determinant V N +\. In other words, the co factor of the element^ in the first row 
is, in fact, the Vandermonde determinant V N +\. Indeed, the cofactor of the nth element in the 
first row, where w = 2, Af+2, is the product of the coefficient a n ,2 in the polynomial 

N 

expansion y = P N (x) = ^a k x k with the Vandermonde determinant V N +\, Furthermore, if x 

k=Q 

and y take on any of the sample values x t and y h for / = 1 , 2, . . . , N 9 iV+l , then two rows of the 
determinant would be the same and the determinant must then vanish. Thus, the requirement 
that the polynomial y = Pjs/(x) pass through the N+l data points (x i9 yj) 9 

N 

y t = P N ( x i) = ^ a k x i >fori= 1,2, TV, N+l, is satisfied. 

k=0 
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For example, a quadratic curve may be found that goes through the sample data set 
(-l,a), (0,6), and (l,c). The three equations are P 2 (-1) = a =a 0 -a\+a 2y P 2 (0) = b=a 0 , and 
P 2 (1) = c =a 0 +ai+a 2 , which imply that b=a 0 , c-a = 2a u and c+a-2b = 2a 2 , so that 



y(x) = P 2 (x) = Z> + x + 



J/ 1 X X 

a 1 -1 1 
Z>100 
c 1 1 1 



c-a c + a-2b 



x 2 , which is also the result of expanding 





1 -1 


1 




a 


-1 


1 




a 


1 


1 




a 


1 


-1 


0 = y 


1 0 


0 


-1 


b 


0 


0 


+ X 


b 


1 


0 


-x 2 


b 


1 


0 




1 1 


1 




c 


1 


1 




c 


1 


1 




c 


1 


1 



the coefficient of 



y being the respective Vandermonde determinant F3 = 2. 



Similarly, a quartic curve may be found that goes through the sample data set (-2 ? a), 

(-1,6), (0,c), (1 5 6), and (2,a). The five equations are P 4 (-2) = a ^ao-2a\+Aa 2 -%a3+\6a^ 

P 4 (-1) = b =ao-ai+a 2 -a3+a 4 , Pa (0) = c =a 0 , P 4 (1) = b =ciQ+a\+a 2 +ai>+a4, and 

P 4 (2) = a =ao+2ai+4a 2 +8a3+l6a4, which imply that c =a 0 , 0 = a\ = a 3 (which also follows 

from the symmetry of the data set), (a-c)-l6(b~c) = -12a 2 , and (a-c)-4(b-c) = 12a 4 , so that 

x a~\6b + \5c 2 a-4b + 3c 4 
y(x) = ^(x) =c x 4- : x 



12 



12 



In various alternative illustrative embodiments, samples may be collected for M data 
points (x/j^/), where i = l, 2, M, and a first degree polynomial (a straight line), 
1 

/>(*) = a 0 +ajX = ^a k x k , may be fit (in a least-squares sense) to the Mdata points (x h yi). 

For example, 100 time data points (M= 100) may be taken relating the degree of sputter 
target consumption, as measured by sputter target life/?, the deposition plasma power/, 
and/or the deposition time T, during one or more metal deposition processing steps, to the 
deposition rate t of metal layers being formed during the respective one or more metal 
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deposition processing steps, resulting in the M data points (p,,/,), (//,/;), and/or (TV,). The 
values may be the actually measured values of the metal deposition processing tool variables 
and/or metal deposition processing parameters, or ratios of actually measured values 
(normalized to respective reference setpoints), or logarithms of such ratios, for example. 
Least-squares fitting is described, for example, in Numerical Methods for Scientists and 
Engineers, by R.W. Hamming, Dover Publications, New York, 1986, at pages 427-443. 

The least-squares criterion may be used in situations where there is much more data 
available than parameters so that exact matching (to within round-off) is out of the question. 
Polynomials are most commonly used in least-squares matching, although any linear family 
of suitable functions may work as well. Suppose some quantity x is being measured by 
making M measurements x„ for /= 1, 2, M 9 and suppose that the measurements x h are 
related to the "true" quantity x by the relation x, =x + s, , for /= 1, 2, M, where the 

residuals 8, are regarded as noise. The principle of least-squares states that the best estimate ^ 
of the true value x is the number that minimizes the sum of the squares of the deviations of 

M M 

the data from their estimate f(4) = = ^X x i •> which is equivalent to the 

j M 

assumption that the average x a , where x a = — V x i f , is the best estimate of the true value x. 
This equivalence may be shown as follows. First, the principle of least-squares leads to the 

M M 

average x a . Regarding / (<g) = ^ ef = ^(*,— £) 2 as a function of the best estimate ^, 



/=i /=i 



minimization with respect to the best estimate £, may proceed by differentiation: 

df( /"^ M MM M 

= -2^ (x, - £ ) = 0 , which implies that £ x,. 4 = 0 x,. -M^ , so that 

d% i= i ,=i ,=i ,=i 
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| M 

= — Vx,. =x fl , or, in other words, that the choice x a = minimizes the sum of the squares 
M 

of the residuals s^. Noting also that — j = ^2 * = > ® » ^ cr i ter i° n f° r a minimum is 

d£> ;=1 

established. 



Conversely, if the average x a is picked as the best choice x a = £, , it can be shown that 
this choice, indeed, minimizes the sum of the squares of the residuals e,-. Set 

M M M M 



/(*J = 5>, -*J =IX -2*„I>, + 2>« -2x a Mx a +Mx 2 a =Y J ^f ~Mx 2 a . If 

/=1 /=1 i=\ i=\ i=\ /=1 

any other value Xb is picked, then, plugging that other value xt> into fix) gives 

M M M M M 

f( x b ) = 2 ( x i ~ x *> ) 2 = 2 x ^ " 2x * 2 *< + 2 x a = 2 x f - 2x t Mx a + • Subtracting f{x a ) 

i=\ /=1 /=1 /=1 /=! 

from gives f(x b ) - ) = m[x] - 2x a x b + x\ ] = M(x a - x b f > 0 , so that X^) 
with equality if, and only if, xt = x a . In other words, the average x a , indeed, minimizes the 
sum of the squares of the residuals e,-. Thus, it has been shown that the principle of 
least-squares and the choice of the average as the best estimate are equivalent. 

There may be other choices besides the least-squares choice. Again, suppose some 
quantity x is being measured by making M measurements x„ for / = 1 , 2, . . . , M, and suppose 
that the measurements x„ are related to the "true" quantity x by the relation x i = x + s , , for 

/ = 1, 2, . . ., M, where the residuals e,- are regarded as noise. An alternative to the least-squares 
choice may be that another estimate % of the true value x is the number that minimizes the 
sum of the absolute values of the deviations of the data from their estimate 

M M 

f (z) - 2KI = 2K ~ 2tl > w hich is equivalent to the assumption that the median or middle 
/=i i=i 
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value x m of the M measurements x h for /= 1, 2, M (if M is even, then average the two 
middle values), is. the other estimate % of the true value x. Suppose that there are an odd 
number M= 2k+l of measurements x h for / = 1, 2, M 9 and choose the median or middle 
value x m as the estimate % of the true value x that minimizes the sum of the absolute values of 
the residuals e,-. Any upward shift in this value x m would increase the k terms \x t -x\ that 

have Xj below x m9 and would decrease the k terms \x i - x\ that have xt above x m , each by the 
same amount. However, the upward shift in this value x m would also increase the term 
\x m - x\ and, thus, increase the sum of the absolute values of all the residuals s,-. Yet another 

choice, instead of minimizing the sum of the squares of the residuals s„ would be to choose to 

x ■ + x 

minimize the maximum deviation, which leads to min max = x midrmige , the midrange 
estimate of the best value. 

Returning to the various alternative illustrative embodiments in which samples may 
be collected for M data points {x h yi), where / = 1 , 2, . . . , M, and a first degree polynomial (a 

straight line), P x (x) = a Q + a x x = ^a k x k , may be fit (in a least-squares sense) to the M data 



points (x,,^/), there are two parameters, ao and a\ 9 and a function F(qo 9 a{) that needs to be 



A=0 



minimized 



as 



follows. 



The 



function 



F(a 09 ai) 



is 



given 



by 



MM M 



a 0 +a l x i —y,f and setting the partial derivatives 




of F(ao,a\) with respect to ao and a\ equal to zero gives 



dF{a 0 ,a x ) 
da 0 




and 



dF(a 0 ,a t ) 
da, 



-2^ i [a 0 +a l x i — xjx, =0, respectively. Simplifying and rearranging gives 
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MM M M M 

a Q M + a x ^ x i y i and a 0 ^ j x i + ^^xf =^x i y i 9 respectively, where there are two 

,=1 ,=l /=! /=] ,=1 



equations for the two unknown parameters a$ and a\, readily yielding a solution. 



As shown in Figure 16, for example, a first degree polynomial (a straight line), 
i 

P x (x) = a Q + a } x = ^a k x k , may be fit (in a least-squares sense) to the M= 5 data points 

k=0 



(1,0), (2,2), (3,2), (4,5), and (5,4). The residuals s h for /= 1, 2, 5, are schematically 
illustrated in Figure 16. The equations for the two parameters ao and a\ are 5ao+15ai = 13 and 
15ao+55ai = 50, respectively, so that, upon multiplying the first equation by 3 and then 
subtracting that from the second equation, thereby eliminating ao, the solution for the 
parameter a\ becomes a\ = 11/10, and this implies that the solution for the parameter ao 
becomes ao = -7/10. The first degree polynomial (the straight line) that provides the best fit, 



in the least-squares sense, is P. (x) = H x = — (- 7 + 1 Ix) , as shown in Figure 1 6. 

H 1 10 10 10 v J 



As shown in Figure 17, for example, a first degree polynomial (a straight line), 



P x (x) = a 0 +a x x = ^a k x k , may be fit (in a least-squares sense) to the M— 1 data points 

A=0 



(-3,4), (-2,4), (-1,2), (0,2), (1,1), (2,0), and (3,0). The residuals e„ for /= 1, 2, 7, are 
schematically illustrated in Figure 17. The equations for the two parameters ao and a\ are 

M M 

a 0 M^-a x Y J x i =7a 0 + a, (-3-2-1 + 0 + 1 + 2 + 3) = (4 + 4 + 2 + 2 + 1 + 0 + 0) and 

/=! /=1 



M MM 

*o£*/ +"i£ x ' 2 =£*i.V,- =«i (9 + 4 + 1 + 0 + 1 + 4 + 9) = (-12 -8-2 + 0 + 1 + 0 + 0), 

/=i /=i /=i 



respectively, which give 7ao =13 and 2Sa\ = -21, respectively. In other words, ao = 13/7 and 
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a\ = -3/4, so that, the first degree polynomial (the straight line) that provides the best fit, in 

13 3 . 
the least-squares sense, is P } (x) — — - — x , as shown in Figure 1 7. 

In various other alternative illustrative embodiments, samples may be collected for M 
data points (x,j>0, where i=l, 2, M, and a polynomial of degree N, 

N 

P N (x) = a 0 + a } x + a 2 x 2 H \-a k x k + — h a N x N = ^a k x k , may be fit (in a least-squares 

sense) to the Mdata points For example, 100 time data points (M= 100) may be taken 

relating the degree of sputter target consumption, as measured by sputter target life p, the 
deposition plasma power/ and/or the deposition time T, during one or more metal deposition 
processing steps, to the deposition rate t of metal layers being formed during the respective 
one or more metal deposition processing steps, resulting in the M data points (/?,//), (fiji), 
and/or (T h ti). The values may be the actually measured values of the metal deposition 
processing tool variables and/or metal deposition processing parameters, or ratios of actually 
measured values (normalized to respective reference setpoints), or logarithms of such ratios, 
for example. In one illustrative embodiment, the degree TV of the polynomial is at least 10 
times smaller than M In various illustrative embodiments, the degree TV of the polynomial 
may be about 2 and the minimum number of data points Mmay be about 25. 

The function F(a^a\^ . .,##) may be minimized as follows. The function F(ao,a\ 9 . . . 5 a^) 

A/ M 

is given by F(a 0 , a x , * • • , a N ) = f = ^ [P N (x i ) - y. ] 2 and setting the partial derivatives of 
F(ao,a\,...,aw) with respect to a 7? for 7 = 0, 1, TV, equal to zero gives 



1=1 /=1 



k=0 



x-j - 0 , for j = 0, 1 , . . . , N, since 
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(xiy is the coefficient of aj in the polynomial P N (x t ) = ^a k xf . Simplifying and rearranging 



gives 



M 



i=\ 



_k=0 



Jt=0 



/=1 



£=0 



M 

\" k S k+J = = Tj , for 7 = 0, 1, N, where 



k=0 



MM N 

^jc* +y = and ^jc/j/,. = T } , respectively. There are equations ^a A 5 ft+> = 7V , for 



7 = 0, 1, also known as the normal equations, for the N+l unknown parameters for 

k=0 9 1, N 9 readily yielding a solution, provided that the determinant of the normal 
equations is not zero. This may be demonstrated by showing that the homogeneous equations 

N 

^a k S k+J = 0 only have the trivial solution ak = 0, for k = 0, 1, ...,]V, which may be shown 
as follows. Multiply the 7th homogeneous equation by aj and sum over all 7, from j = 0 to 

N N N N M M f N \f N \ M 

7=0 *=0 y=0 *=0 i=l /=1 V*=0 J\J=0 J i=\ 

which would imply that P/v(x,) = 0 ? and, hence, that ak = 0 9 for k=0 9 1, . .., yV, the trivial 
solution. Therefore, the determinant of the normal equations is not zero, and the normal 
equations may be solved for the 7V+1 parameters a*, for k = 0, 1 , ...,N 9 the coefficients of the 



least- squares polynomial of degree N, P N (x) = ^a k x k , that may be fit to the Mdata points 



N 



Finding the least- squares polynomial of degree N 9 P N (x) = ^& k x k , that may be fit to 

k=0 

the Mdata points (x h yi) may not be easy when the degree Nof the least-squares polynomial is 

N 

very large. The N+l normal equations ^a k S k+J - T jy for 7 = 0, 1, N 9 for the JV+l 
unknown parameters a k9 for k = 0, 1, ...,JV } may not be easy to solve, for example, when the 
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degree N of the least-squares polynomial is much greater than about 10. This may be 
demonstrated as follows. Suppose that the M data points {x h y t ) are more or less uniformly 

distributed in the interval 0 < x < 1 , so that S k+i = Yx* +J « M [x k+J dx = . The 

k+J U h k + j + l 

resulting determinant for the normal equations is then approximately given by 



J k+j\ 



M 



k + j + 1 



M 



N+\ 



1 



* + j + 1 



= M N+1 H N+] ,forj,k=0, 1, ...,N, where H N ,forj 9 k=0, 1, 



N-\, is the Hilbert determinant of order N, which has the value 
[0!l!2!3!--(iV-l)!] 3 



N\(N + l).(N + 2).---(2N - 1) 



that approaches zero very rapidly. For example, 



1! ll' 2 2!3! 



1 * 



2 3 



= i_i = j_ jand „ j= M* 

3 4 12 3 3!4!5! 



1 1 1 

1 2 3 

111 

2 3 4 

111 

3 4 5 



5 where 



1 1 1 

1 2 3 

111 

2 3 4 

111 

3 4 5 



2 3 


1 


1 2 


1 

+ - 


1 2 


1 1 

3 4 


4 


1 1 

3 4 


5 


1 1 

2 3 



-I_L_IJ_ + IJ_- 1 1 



3 72 4 12 5 12 216 240 2160 



. This suggests 



that the system of normal equations is ill-conditioned and, hence, difficult to solve when the 
degree TV of the least-squares polynomial is very large. Sets of orthogonal polynomials tend to 
be better behaved. 



As shown in Figure 18, for example, a second degree polynomial (a quadratic), 



P 2 (x) = a Q -\-a x x + a 2 x 2 =^a k x k , may be fit (in a least-squares sense) to the M=7 data 

k=0 

points (-3,4), (-2,2), (-1,3), (0,0), (1,-1), (2,-2), and (3,-5). The residuals e h for /= 1, 2, 7, 

2 

are schematically illustrated in Figure 18. The three normal equations are ^a k S k+J = T J9 for 



A=0 



7 = 0, 1,2, where = S k+J and ^xfy f s T } , respectively, for the three parameters ao, 



i=i 
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and a\ and a 2 . This gives ^a k S k =T 0 , ^a k S k+} =T l9 and ^^^ +2 = T 2 5 where 

k=0 k=0 k=0 

S 0 = J>, 0 =7;S X = £x, =(-3-2-l + l + 2 + 3)=0;S 2 = ^> 2 =(9 + 4 + 1 + 1 + 4 + 9)- 28; 
/=i /=i /=i 

S 3 = ^>, 3 =(-27-8-l + 0 + l + 8 + 27) = 0;S 4 = ^>, 4 =(81 + 16 + l + 0 + l + 16 + 8l) = 196; 
/=i /=i 

T 0 =Y j y i =(4 + 2 + 3 + 0-l-2-5) = l;7; = J>,.j>/ =(- 12 - 4 - 3 + °- 1 - 4 - 15 )^ ~ 39 ^ 
/=i /=i 

7 

and T 2 = ^ jc 2 ^,. =(36 + 8 + 3 + 0-1-8 - 45) = -7 , so that the normal equations become 

2 2 
J^S* = T 0 = 1 = 7 «0 + 0a i + 28 ^2 = 7a 0 + 2 8^2 > = ^ = " 39 = ^0 + ^1 + 0fl 2> 



A=0 



and ^a k S k+2 = T 2 = -7 = 28a 0 + 0a } + 196a 2 = 28a 0 + 196a 2 , respectively, which imply 



A=0 



(upon multiplying the first normal equation by 7 and then subtracting that from the third 
normal equation) that -14 = -21ao ? that 28ai =-39 (from the second normal equation), and 
(upon multiplying the first normal equation by 4 and then subtracting that from the third 
normal equation) that -1 1 = $4a 2 , giving 3a 0 = 2, 28a { = -39, and 84a 2 = -1 1, respectively. In 
other words, a 0 = 2/3, a\ = -39/28, and a 2 = -1 1/84, so that, the second degree polynomial (the 
quadratic) that provides the best fit, in the least-squares sense, is 

p (x) - — — — x- — x 2 = — (56-117x-llx 2 ), as shown in Figure 18. 
2W 3 28 84 84 V 7 

As shown in Figure 19, for example, a second degree polynomial (a quadratic), 

2 

P 2 (x) = a 0 + a x x + a 2 x 2 = ^a k x k , may be fit (in a least-squares sense) to the M— 6 data 

*=0 

points (0,4), (1,7), (2,10), (3,13), (4,16), and (5,19). The residuals e„ for /= 1, 2, 6, are 
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schematically illustrated in Figure 19. The three normal equations are ^a k S k+J =T J9 for 

k=0 

6 6 

j = 0, 1,2, where ^xf +J = S k+J and ^ */>>,■ = 7,. , respectively, for the three parameters a 0 > 
/=! /=i 

2 2 2 

and 0i and #2. This gives ^a k S h =T Q9 ^a k S k+l =T l9 and ^a k S k+2 =T 2 , where 

k=0 k=0 k=0 

S 0 = =6:5, = =(0 + l + 2 + 3 + 4 + 5) = 15;S 2 = =0 + 4 + 9 + 16 + 25)- 55 ; 

S 3 =^xf =(0 + 1 + 8 + 27 + 64 + 125)= 225;S 4 = =(0 + 1 +16 + 81 + 256 + 625) = 979; 

#=i /=i 

7 0 = Y 4 y i = (4 + 7 + 10 + 13 + 16 + 19)= 69;7; = Y d x i y i =(° + 7 + 20 + 39 + 64 + 95 ) = 225 ^ 

/=1 /=! 
6 

and J 2 = ^ x, 2 ^,- =(0 + 7 + 40 + 1 17 + 256 + 475) = 895 , so that the normal equations become 

2 2 

Y, a k s k = T o =69 = 6a 0 +15a, + 55a 2 , = r i = 225 = l5a ° + 55a ' + 225a 2> 



k=0 



^ a k S k + i = T i = 895 = 55a 0 + 225^ + 979a 2 , respectively, which imply (upon multiplying 



the second normal equation by 4 and then subtracting that from the first normal equation 
multiplied by 10) that -210 = -70#i-350a2 ? and (upon multiplying the second normal equation 
by 11 and then subtracting that from the third normal equation multiplied by 3) that 
210 = 70ai+66a2. However, adding these last two results together shows that 0 = ^2. 
Furthermore, 3 = ct\. Therefore, using the fact that 3 = a\ and 0 = the normal equations 

2 2 

become JX 5 * =T o =69 = 6a 0 +45, JX S a+i = t i =225 = 15o 0 +165, and 

*=0 k=0 



^a k S k+2 = T 2 = 895 = 55a 0 + 675 , respectively, which all imply that 4 = a 0 - In other words, 

k=0 

ao = 4, a\ = 3, and <?2 = 0, so that, the second degree polynomial (the quadratic) that provides 
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the best fit, in the least-squares sense, is P 2 (x) = 4 + 3x + Ox 2 = 4 + 3x , which is really just a 
straight line, as shown in Figure 19. The residuals e„ for z = 1, 2, 6, all vanish identically 
in this case, as schematically illustrated in Figure 19. 



In various other alternative illustrative embodiments, samples may be collected for M 
5 data points (x i: yi), where i= 1,2, . . ., M, and a linearly independent set of 7V+1 functions j£(x), 

N 

fory = 0, 1, 2, N, Xx) = flJ 0 W + ^ 1 W + --- + a//i) + - + a^(x) = ^a/ j (x), 

7=0 

may be fit (in a non-polynomial least-squares sense) to the Mdata points (x^yi). For example, 
100 time data points (M= 100) may be taken relating the degree of sputter target 
'if consumption, as measured by sputter target life p, the deposition plasma power/ and/or the 

jj si 

£V 10 deposition time 7, during one or more metal deposition processing steps, to the deposition 

p rate / s of metal layers being formed during the respective one or more metal deposition 

y i 

a processing steps, resulting in the M data points (p i9 ti) 9 (fiji), and/or (TV,)- The values may be 

Q 

In the actually measured values of the metal deposition processing tool variables and/or metal 

2f deposition processing parameters, or ratios of actually measured values (normalized to 

15 respective reference setpoints), or logarithms of such ratios, for example. In one illustrative 
embodiment, the number N+l of the linearly independent set of basis functions fj(x) is at 
least 10 times smaller than M. In various illustrative embodiments, the degree N of the 
polynomial may be about 2 and the minimum number of data points Mmay be about 25. 



The function F(a 0 ,ai,.. .,a N ) may be minimized as follows. The function F{a^a\^..,a^l) 

M M 

20 is given by F(a Q ,a x ,- -,a N ) = ^£, 2 = ^ [>>(*,■)- >V P and setting the partial derivatives of 
F(ao,a\,...,a N ) with respect to a J? for 7 = 0, 1, N 9 equal to zero gives 
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k=0 



f J (x i ) = 0 9 fOTj = 0 9 1, 



since fj(x?) is the coefficient of aj in the representation j>(x f .) -^a k f k (x,) . Simplifying gives 



/=1 



.Jt=0 



/j(*,")=Z fl * 



*=0 



M 

£/*(*,)//*,) 



/=1 



*=0 



S 2>A,, = - s7 *y. for 7 = 0, 1, 

k=0 ;=1 



N, where 2 AC*/)//*/) = 5 m Z/ y (>, ).y, = ^ 7 , respectively. There are N+l 

N 

equations Z ci k S k } = Tj , for j = 0 9 1 , ...,iV 5 also known as the normal equations, for the 7V+1 



unknown parameters a*, for k = 0, 1 , . . . , TV, readily yielding a solution, provided that the 
determinant of the normal equations is not zero. This may be demonstrated by showing that 

N 

the homogeneous equations ^ J d k S kJ = 0 only have the trivial solution ^ = 0, for k — 0, 1, 

k=0 

... 9 N 9 which may be shown as follows. Multiply the yth homogeneous equation by aj and sum 

N N N N M M f N \f N 

over ally, Y, a jH a k S kj = Z a ; Z a * Z A ( x / >// x < ) = Z Z a */*( x /) Z a j// X <) 

^ M 



V N 



VJ=0 



, but 



M f N Y 1 A-? 

Z Z a *A(^/) Z a >/y( x ') =Z0 ; ( x /)) 2=0 ' which would ^P 1 ^ that ^( x <) = °> and > 

i=l \k=0 \j=o J i=\ 



hence, that a* = 0, for k = 0, 1 , . . . , N 9 the trivial solution. Therefore, the determinant of the 
normal equations is not zero, and the normal equations may be solved for the N+l parameters 
a k , for k = 0 9 1, N 9 the coefficients of the non-polynomial least-squares representation 

N 

y( x ) ~ ^, a jfj( x ) ^at may be fit to the Mdata points (x h yi) 9 using the linearly independent 



7=0 



set of N+l functions fj(x) as the basis for the non-polynomial least-squares representation 



jK*) = Z a ;/>(*)- 
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If the data points (x h yi) are not equally reliable for all M 9 it may be desirable to weight 
the data by using non-negative weighting factors w,-. The function F(a 0 ,a u ...,a N ) may be 
minimized as follows. The function F{ao 9 a\ 9 ,., 9 a^) is given by 



IVl IVl 

F(a Q 9 a x , • • • , a N ) = ^ w,.e f — w,. [y (x, ) - >>,- ] 2 and setting the partial derivatives of 

,-=l ;=i 

F(a 0 ,a\,...,a N ) with respect to a,-, for J = 0 9 1, equal to zero gives 



da, 



i=\ 



/,-(*,) = 0, for y = 0, 1, 



N, since fjipci) is the coefficient of ctj in the representation j(x f . ) = ^ a k f k (x, ) . 



jfc=0 



Simplifying gives 



*=0 



A=0 



^Y J a kS k j^Y J ^ i f J {x i )y i ^T j for y = 0, 1, TV, where 



^ w //a (*/)/,■(*,•) — $kj ^ WifjiXf))?,- = Tj , respectively. There are 7V+1 equations 

/=i /=i 

N 

^a a k^kj —T J9 for 7 = 0, 1, N 9 also known as the normal equations, including the 



A=0 



non-negative weighting factors w i9 for the 7V+1 unknown parameters for k~ 0, 1, ... 9 N 9 
readily yielding a solution, provided that the determinant of the normal equations is not zero. 

N 

This may be demonstrated by showing that the homogeneous equations ^o k S kJ =0 only 

A=0 

have the trivial solution a* = 0, for k = 0, 1 , . . N 9 which may be shown as follows. Multiply 
the yth homogeneous equation by a, and sum over all j 9 



N N N N M M f N \f N ^ 

H a jH a k s *j = ll a jll a kYa w ifki x i)fji x ^ = Y, w ' H a kfki x i) Y. a jfA x i"> 

>=0 Jk=0 7=0 A=0 /=! /=! V*=0 J\ J=0 



but 
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M ( N \f N \ M 

Z w / Z a kfk o, ) Z °y/y <x ) = Z w i (y(*i )) 2 = 0 > which would ^p 1 ^ that y s °> ^ 



V*=o 



hence, that a* = 0, for £ = 0 5 1 5 ... 9 N, the trivial solution. Therefore, the determinant of the 
normal equations is not zero, and the normal equations, including the non-negative weighting 
factors w h may be solved for the N+l parameters for k = 0, 1, . . . , TV, the coefficients of the 

N 

5 non-polynomial least-squares representation y(x) = Z a y//00 that may be fit to the Mdata 

7=0 

points (xiji), using the linearly independent set of N+l functions fj(x) as the basis for the 

N 

non-polynomial least-squares representation = Z a y/, ( x ) ? including the 

yp non-negative weighting factors Wj. 

H Any of the above-disclosed embodiments of a method according to the present 

O 10 invention enables the use of parametric measurements sent from measuring tools to make 

y s 

JL supervisory processing adjustments, either manually and/or automatically, to improve and/or 

l fJ better control the yield. Additionally, any of the above-disclosed embodiments of a method of 

p, manufacturing according to the present invention enables semiconductor device fabrication 

with increased device accuracy and precision, increased efficiency and increased device 
15 yield, enabling a streamlined and simplified process flow, thereby decreasing the complexity 
and lowering the costs of the manufacturing process and increasing throughput. 

The particular embodiments disclosed above are illustrative only, as the invention 
may be modified and practiced in different but equivalent manners apparent to those skilled 
in the art having the benefit of the teachings herein. Furthermore, no limitations are intended 
20 to the details of construction or design herein shown, other than as described in the claims 
below. It is therefore evident that the particular embodiments disclosed above may be altered 
or modified and all such variations are considered within the scope and spirit of the invention. 
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In particular, every range of values (of the form, "from about a to about 6," or, equivalently, 
"from approximately a to 6," or, equivalently, "from approximately a-b") disclosed herein is 
to be understood as referring to the power set (the set of all subsets) of the respective range 
of values, in the sense of Georg Cantor. Accordingly, the protection sought herein is as set 
forth in the claims below. 
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